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The expectation value of a Wilson loop in a Chern-Simons theory is a knot invariant. Its 
skein relations have been derived in a variety of ways, including variational methods in which small 
' deformations of the loop are made and the changes evaluated. The latter method only allowed to 

0^ ' obtain approximate expressions for the skein relations. We present a generalization of this idea 

G^ : that allows to compute the exact form of the skein relations. Moreover, it requires to generalize the 

' ' resulting knot invariants to intersecting knots and links in a manner consistent with the Mandelstam 

O , identities satisfied by the Wilson loops. This allows for the first time to derive the full expression 

D ■ for knot invariants that are suitable candidates for quantum states of gravity (and supergravity) in 

the loop representation. The new approach leads to several new insights in intersecting knot theory, 
QQ , in particular the role of non-planar intersections and intersections with kinks. 

(N ; 
> ■ 

^ ' I. INTRODUCTION 

^ : 

. Witten realized some years ago that the expectation value of a Wilson loop W{'^) in a Chern-Simons theory 
' was a knot invariant. This follows from the fact that Chern-Simons theories are diffeomorphism invariant and that 
, the Wilson loops are observables for such theories, having therefore diffeomorphism invariant expectation values. The 
' resulting knot invariant is the Kauffman bracket for the case of an SU (2) Chern-Simons theory and for the case 
^ ' of SU{N) is a regular isotopic polynomial associated with the HOMFLY H polynomial. These results were derived 
• based on the calculations of Moore and Seiberg |4| for the monodromies of rational conformal field theories. The 
O i' knowledge of the Yang-Baxter relations satisfied by the monodromies translates immediately into skein relations for 
D the polynomial in question. 

Independently, Smolin and later Cotta-Ramusino, Guadagnini, Martellini and Mintchev |^ noted that a simpler 
J> , heuristic derivation of the skein relations was possible. The idea is similar to the Makeenko-Migdal approach 
to Yang-Mills theories. It is based on studying the changes in the expectation value of the Wilson loop when one 
. , performs small deformations. This calculation can be done explicitly to first order in the deformation. The results 
■ can be interpreted as skein relations to first order in the inverse coupling constant of the theory, which is tantamount 
■ ■ ' to determining the knot polynomial to first order. 

The latter method is quick and computationally efficient, and has a simple generalization to the case of intersecting 
loops 1^,^. The main drawback, especially in the case where the result is not known by other methods, is that one 
only gets the skein relations to first order in the inverse of the coupling constant of the theory. It is therefore of 
interest to find a suitable generalization that would yield the exact skein relations to all orders. This is the main 
purpose of this paper. 

On the other hand, the subject of intersecting knot invariants has received little attention and is of paramount 
importance for the construction of quantum states of gravity in the loop representatio n Fo| JTl[|. In this approach, 
based on the canonical quantization of general relativity in terms of Ashtekar variables |1 jj^wavefunctions are knot 
invariants due to the diffeomorphism symmetry of general relativity Jl3t . The Hamiltonian constraint has only a 
non-trivial action at intersections jlj]. Only intersecting knots are associated with non-degenerate spacetimes [ [l5| . 
Whenever one generalizes an invariant of smooth loops to take values on intersecting loops there is generically freedom 
in how the invariant is defined, as long as it is compatible with the Reidemeister moves. However, in the case of 
quantum gravity, wavefunctions have to be compatible with a set of constraints among functions of loops known 
as the Mandelstam identities. These identities naturally involve intersecting loops and severely limit the possible 
generalizations of invariants to intersections. We will show in this paper how to generalize the invariants stemming 
from Chern-Simons theory to be compatible in an exact way with the Mandelstam identities. This in particular also 
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defines the values of the invariants for multicomponent links. This is of particular relevance for quantum gravity 
since it is known that the exponential of the Chern-Simons form built from the Ashtekar connection is an exact 
solution to all the constraints of quantum gravity If one wishes to find the counterpart of this state in the 

loop representation one ends up computing exactly the same integral as the expectation value of a Wilson loop in a 
Chern-Simons theory. 

One additional motivation for the construction we present is that the Chern-Simons state not only arises in canon- 
ically quantized vacuum general relativity but also in other contexts, like Einstein- Yang-Mills theories |16| and 
supergravity [l^ . In these cases the (super)gauge group of the associated Chern-Simons theory differs from that of 
gravity and therefore so do the resulting invariants. In the particular case of supergravity the resulting invariant had 
not been computed by other means, and turns out to be associated with the Dubrovnik-Kauffman polynomi al |19|| . 

A first attempt to obtaining a finite prescription from variational calculations was made by Briigmann pG]. In 
particular, the idea of exponentiating the infinitesimal transformation was presented there. Because of ambiguities 
in the formulation presented in that paper, one could only check, given the exact skein relations, that they were 
compatible with the formulation. Here we add two key elements that make the construction a well defined prescription: 
on the one hand we offer a justification of why can the infinitesimal results be exponentiated; on the other hand we 
make crucial use of the Mandelstam identities to uniquely fix a prescription for the exponentiation. The prescription 
now allows, in a case where one does not know the result beforehand, to compute the resulting polynomial. 

The organization of this article is as follows. In the next section we will discuss, using the non-Abelian Stokes 
theorem, how to perform a finite deformation of the expectation value of the Wilson loop. In section III we discuss the 
deformation of twists and kinks and in section IV of planar and non-planar intersections. We discuss the implications 
of the results in section V. 



II. THE NON-ABELIAN STOKES THEOREM AND FINITE DEFORMATIONS OF LOOPS 



One is interested in establishing skein relations for the following function of a loop, 

< W{-i) >= J dAexp (^^Scs^ W{j) 

where 



IF (7) = Tr 



Pexp [i jdy" Aa{y) 



(1) 



(2) 



where 7 is a closed curve in a three manifold, Aa is a connection in a semi-simple Lie algebra and Scs is the action 
of a Chern-Simons theory. 



Scs — 



d^xe^^^Tri^.afc^, - -lAaA^A,] 



(3) 



Finding skein relations involves relating the value of the < W{'^) > for different loops. These loops are determined 
by the replacement of over-crossings by under-crossings in a planar projection. Skein relations also arise between a 
loop with and without a twist, over-crossings and intersections. The usual notation for the several types of crossings 
is shown in figure |l|. We will see later on that the consideration of invariants with intersections requires several 
other crossings. For instance, the skein relations that define the Kauffman bracket knot polynomial on loops without 
intersections are. 



KiL+) 
i^(unknot) — 1. 



(4) 
(5) 
(6) 
(7) 



With these relations the polynomial is completely characterized for any link. What Witten showed using 
conformal field theory techniques is that i < 11^(7) > satisfies the above relations. Here we will show it by performing 
directly deformations of the loops on the expression of the expectation value. 

In order to do this, we need to study deformations of Wilson loops. Wilson loops are traces of holonomies. It turns 
out that the information needed to compute a holonomy is less than that present in a closed curve. Several closed 
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FIG. 1. The different crossings involved in the skein relations for invariants of non-intersecting loops 

curves yield the same holonomy. In this paper we will use the word "loop" to denote the equivalence class of curves 
that yield the same holonomy for all connections |^. Loops form a group structure called the group of loops . 

It is well known that if one adds to a loop another loop of infinitesimal area, the change in the Wilson loop can 
be coded in terms of an infinitesimal operator in loop space called the loop derivative . This is shown in figure |^. 
The concrete expression is, 



o 57) = (I + -a'^''A,,«))T4^(7) 



(8) 



where o denotes composition of loops, cr°^ is the infinitesimal element of area of the the loop 5^ and tt^ is an open 
path connecting the basepoint of the loop 7 to the point x at which one adds the infinitesimal loop with infinitesimal 
element of area cr"''. Notice that the loop derivative Aah depends on the path used to compute it and we denote so 
in its expression. The loop derivative is related to the infinitesimal generators of the group of loops |I6|. 




FIG. 2. The infinitesimal loop that defines the loop derivative. 

One can write [ p^ an expression for an operator U{"f) that adds a finite loop 7 to a Wilson loop in terms of the 
loop derivative. Let 7(5) be a parameterized curve belonging to the equivalence class defining the finite loop 7 with 
s G [0, 1]. Consider a one-parameter family of parameterized loops 77(5, t) interpolating smoothly between 7(3) and the 
identity loop, such that 77(3, 0) is in the equivalence class of the identity loop and rj{s, 1) = 7(s). Consider the curves 
77(5, 1) (= 7(5)) and 77(5, 1 — e). The two curves are drawn in figure || and differ by an infinitesimal element of area. 
The whole purpose of our construction will be to cover the infinitesimal area separating the two mentioned curves 



^ Other authors call these objects "hoops" to denote "holonomic loops" |21 
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FIG. 3. The construction of a finite loop from the loop derivative. The curves 5rji are determined by two elements in the 
family rj{t). 

with a "checkerboard" of infinitesimal closed curves such that along each of them one can define a loop derivative. 
One can therefore express the curve 7(5) as 

7(5) = lim 77(5, 1 — e) o 5r]i o • • • o 5r]n, (9) 

n — >QC 

where the Sr/i are shown in figure |^. Analytically, in terms of differential operators on functions of loops we can write^ 



(10) 

where 'f]{t,s) = dr]{t, s)/dt and ri'(t,s) = drjit, s) / ds. It is immediate to proceed from rjil — e,s) inwards just by 
repeating the same construction, and so continuing until the final curve is the identity. The end result is 



vI/(7) = Texp( 1^ dtj^ dsri%t,s)rj'\t,s)Kb{ri{t)l)]^{rim = U{l)^{rim 



(11) 



where the outer integral is ordered in t (T-ordered) . This result is the loop version of the non- Abelian Stokes theorem 
of gauge theories and it shows that the loop derivative is a generator of loop space, i.e., it allows us to generate any 
finite loop homotopic to the identity. Due to the properties of the group of loops |16|, the construction is independent 
of the particular family of loops used to go from the identity element to the final loop 7. 
It is useful to rewrite the expression of the operator U as. 



where 



C/(7) = Texp / dt5{t) 
Jo 



(12) 



(13) 



The operator S is closely related to the unparameterized "connection derivative" |16[| . 

We would now like to apply the above deformation to a portion of a loop in order to construct the elements that 
appear in the skein relations. 



We drop the s dependence of 77 where it is not relevant. 
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FIG. 4. Adding a twist to a loop using the finite deformation operator 



III. SKEIN RELATIONS ASSOCIATED TO TWISTS AND KINKS 



A. Twists 



We start by the simplest skein relations, those that relate the value of the invariant with and without a twist. 
Starting from a regular portion of a loop 77^ going from the origin to x, one can add a twist by considering a family 
of loops of the form, 



(14) 



where u'^{s) is a vector along the loop that materializes the deformation shown in figure ^ 

To compute the deformation, wc first evaluate the action of the loop derivative on the expectation value of the 
Wilson loop. 



(15) 



where we have taken into account that the deformation is along one element of the the family of loops ri{t). are 
the generators of the algebra su{2). In this calculation wc have assumed that the loop derivative acts on < W{'^) > 
by simply acting inside the functional integral ^^(7). This is a strong hypothesis, notice that the integral < 1^(7) > 
is diffeomorphism invariant and therefore the limit involved in the definition of the loop derivative is singular. We are 
assigning a value to that limit by permuting the functional integral and the limit involved in the loop derivative. 

We now use the fact that the field tensor is the dual of the magnetic field, which can in turn be obtained through 
a functional derivative of the exponential of the Chern-Simons state, 



\Scs 



(16) 



(17) 



Substituting this expression in (|15|) and integrating the functional derivative by parts we get, 

and acting with the functional derivative on the holonomy we get, 
a^'AMtTo) < W{j) >=-^J dAo'^^e^'^^^tabc dy^S^x - 2/)Tr (r'ff (,7(t)°)iJ(7„^)r'ff (7°)iJ(r;(i)^)) . (18) 

Notice that the previous expression is distributional, involving a one-dimensional integral of a three-dimensional 
Dirac delta. It is remarkable that one can use it to obtain an expression for the operator U that generates finite 
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deformations. In order to see this, we first write an expression for the operator 5 associated with the deformation 
induced by the vector field u(s), 



Suit) < 14^(7) dsu^is)iv'\s)+tu'\s))eabc <f dy^6Hvis,t) ~ y) {Tt {r' H{T^{i)l)H{^yy H{^l)H{r^{t)l))) 

(19) 

In order to construct the finite deformation operator we need to exponentiate 5. As is shown in appendix A, one 
can find a regularization such that 5's at different points commute when acting on the expectation value of the Wilson 
loop in Chern-Simons theory (in general they do not commute, see ]l6[|). Therefore the T-ordered exponential reduces 
to an ordinary exponential and we get, 

t/(7) = exp(^^ 5z{t)dt^ (20) 
when acting on < W{'^) > and the integral in the exponent can be computed explicitly, 

dt6s{t) < W{-/) T— < Wi^) > (21) 

and the sign =F corresponds to the kind of crossing generated in the loop by the vector field u. Positive sign corresponds 
to the right hand rule. To go from ^ to the last expression we have made use of the Fierz identity for SU{2), 

kA kC _ IxA rC _ Ix^xC 

and we performed explicitly the three one-dimensional integrals. The result is the normalized oriented volume sub- 
tended by the deformation. Strictly speaking, this result implies a choice of regularization, since the type of integral 
that one is left with is of the form dtS{t) (see for details). This choice in the regularization is tantamount to 
introducing a new parameter in the derivation, corresponding to the value of the above integral. Throughout this 
paper we will take it to be unity. Otherwise, it would imply a multiplicative shift in the value of k. 

The above expression can be summarized in terms of the notation for skein relations we introduced before as, 

< W{L+) > = cxp f - — j < W{Lo) > (23) 

< W{L^) > = exp (^) < WiD > (24) 



which coincide with the skein relations for the Kauffman bracket (^||5|) if one makes the identification 

(Z = exp(- — ). (25) 

It is well known that perturbative techniques like the ones we are using here fail to capture the additive shift in the 
coupling constant first observed by Witten [Q. A discussion, and a proposal to amend perturbative techniques to 
capture this effect through a recourse to the semi-classical approximation can be found in Awada ||2^. We could 
proceed in the same fashion here and modify the value of the coupling constant, but we will leave it as it is for 
simplicity. 



B. Kinks 



A kink (discontinuity of the tangent) is a diffeomorphism invariant feature of a loop. We will now show that the 
expectation value of a Wilson loop in a Chern-Simons theory is not sensitive to the presence of kinks in the loop 
under the kind of regularization we are using in our calculations. This is an important result because it directly relates 
to the value of the invariant when one has intersections with kinks, which are crucial to implement the Mandelstam 
identities. We will discuss this at the end of this section and will see later that the choice we make for the treatment 
of the kinks is central to establishing consistency. 
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FIG. 5. Deforming a kink into a smooth section 
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FIG. 6. Double intersections without and with kinks 

Let us consider a loop with a kink like that shown in figure |^. The kink defines a plane. We can deform the kink 
into a smooth section of the loop through a deformation in the plane. The calculation is exactly the same as that of 
a smooth section, so we will not repeat it here. The only precaution is to consider the appropriate tangent vector 
at both sides of the kink. The reason why we do not need more details is that since the deformation is planar, the 
contribution vanishes. Basically one gets a contribution similar to ([l9| ) that has three coplanar vectors contracted 
with the Levi-Civita symbol. 

The case of intersections with kinks is essentially similar. There are two different types of (double) planar intersec- 
tions with kinks, as shown in figure ^. 

In all cases the operator U{j) deforms independently each line at the intersection. For the cases with kinks that 
means that each kink can be deformed into a smooth section. That means that in cases L\y and L^^ the intersection 
can be removed by the deformation and we get the skein relations 



< W{Lw) > = < W{Lo) > 



(26) 
(27) 



that are shown in figure 

Notice that in all cases the removal of kinks at planar intersections can be accomplished given a straightforward 
regularization of the deformation operator. This provides a justification for the identification of Lyy with Lq that has 
been used in other works |2^,|2^. We will return to discuss intersections with kinks at the end of the next section, 
where we will analyze a deformation that produces double lines in the loop that is of interest in the context of the 
exponentiation of the skein relation for straight through intersections, which we discuss now. 



IV. SKEIN RELATIONS FOR INTERSECTIONS 



FIG. 7. Skein relations for intersections with kinks 
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FIG. 8. Deformation of a "straight through" intersection. The deformation vector u is perpendicular to the plane determined 
by the intersection 



A. Infinitesimal skein relations 



In the case of a "straight through" intersection (no kinks), the calculation is different than in the previously discussed 
cases. In this case, the Wilson loop is the trace of the product of the holonomies along the petals defined by the 
intersection. For instance in the case of figure ^ one could write it as 1^(7) — Tr(iJi2i?34). The action of the loop 
derivative on the holonomy is independent of the presence of the intersection, the expression is exactly the same as 
(p^). The difference in the construction arises when one integrates by parts. Since the loop has a double point at the 
intersection, the functional derivative with respect to Aa has two contributions at that point, corresponding to its 
action on the holonomy when it traverses that point the first and the second time. One of the contributions vanishes 
since it produces a term proportional to the tangent of the loop in the same plane as the deformation and therefore 
spans no volume. The other contribution is the one of interest. It can be written as, 

a^^'AMtYo) < W{^) > = [ dAa'^^'e^^^^abc I dy^5\x - y) (28) 

xTr (T'i7(ry(t)°)i/34(7o')^ri2(7Dr'i?i2(7,0^f34(7?)^^(^(0S)) ■ 

where we have assumed that the origin o of the loop is in the petal 34. We now use the Fierz identity ( p2| ) and 
compute the integral of the 5 operator, 

dt5n{t) < Tr(i/i2(7)i^34(7)) >= -- ^Tr(i/i2(7)i?34(7)) + — Tr(Hi2(7))Tr(i/34(7)) (29) 

Which can be rewritten as 

^ 2?7r 
dtSsit) < W{Li) -— < W{Li) > +— < W{Lw) > ■ (30) 

A couple of remarks are in order. First, there is a difference between equations (^) and (^o|), since equation (|29| ) 
refers to the whole loop and equation ( |3C| ) only to its intersection. The former has information about the connectivity 
of the loop, the latter does not. Therefore, if one is to claim that equation (^) follows from equation ( p9|) one 
should prove that it does so for any connectivity. For the case of double intersections, there are only two different 
possible connectivities. It is straightforward to check that equation (|30| ) holds independently of the connectivity 
chosen. Second, in spite of similarities to equation (|^), equation ( |30| ) is only valid to first order in 1/fc and needs to be 
exponentiated to obtain the skein relation. Since equation ( ^0| ) shows that the action of S mixes L/ and Lw, in order 
to exponentiate it we need to compute the action of the deformation operator on Lw- Notice that this is a different 
deformation than the one computed in the last section, which was coplanar with the intersection. We discuss of the 
resulting calculation at the end of this section. 



Before doing that computation, it is worthwhile observing that one could combine the Fierz identity (22) with the 
following identity, 

e^^esD^Sp%^6^S^ (31) 

to get, 

^ D ^ J\.^d(>b - ^ ^bd): (32) 
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FIG. 9. The skein relation for intersections requires the analysis of this type of deformation of an intersection with kinks. 



and using this expression in (|2^) we get 

rl 

dtSnit) < Tr(i/i2(7)iJ34(7)) >= yTr(Fi2(7))T>(i?34(7)) + y Tr (7712(7)^3"/ (t))- (33) 

In contrast to (|2^) this expression leads to different first-order skein relation depending on the connectivity of the 
knot. If the the connectivity is such that the original loop has a single component, then the skein relation is, 

/ tTT ITT 

J dtSsit) < W{Lj) >= y < W{L„) > +y < W{Lo) > . (34) 

where the element Lq) is defined in figure |l|. The interest of this skein relation stems from the fact that the original 
definition of the Kauffman bracket was given in terms of relations of this kind. Moreover, Major and Smolin |2j] 
used this identity to derive the binor identity from the Kauffman bracket. The original definition of the bracket differs 
from the one we use here in a factor (—1) to elevated to the number of connected components of the loop, which fixes 
the connectivity difference we encountered above. 

To conclude, we discuss the deformation of an Lw as mentioned before. We need to compute the deformation of 
an intersection with a kink that is shown in figure |^. The loop starts at the origin, goes through the deformation 
determined by the vector u, which we denote with a dotted line and then traverses the petal 12 going through the 
kink and ends with the petal 34. The resulting loop has double lines. In ordinary knot theory double lines are not 
considered. They could be incorporated through additional sets of skein relations, as we discuss in subsection F. If one 
pursues a calculation similar to the ones we have been doing up to now for this case, one encounters two contributions, 
stemming from the volumes spanned by the deformation and the tangents to the loop at the intersection in lines 2 or 3. 
The integrals involve terms of the form 5^{y)Q{y). The Heaviside function Q{if) limits the integral in the contribution 
to the corresponding petal. Although these kinds of expressions can be regularized, there does not appear any natural 
way of assigning relative weights to the two contributions from the petals. The way we handle it here is to write the 
contribution in terms of two arbitrary factors a and (3, 

ITT 

/ dt<W{Lw±)>=±^{a<W{Li)>+[3<W{Lw)>). (35) 
Jo ^ 

where we have denoted the elements that result from doing the deformation up or down with respect to the plane 
determined by the intersection. 

We will see in the next section that the arbitrary factors can be uniquely determined by the Mandelstam identities 
when we exponentiate the operator. 

B. Exponentiation and Mandelstam identities 

The skein relations involving L± and Lj are obtained studying the deformation of an intersection. In order to 
obtain these deformations in a finite form we need to exponentiate the differential expressions we obtained in the 
previous section. Specifically, we have 
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FIG. 10. The different loops that arise in the Mandelstam identities involved in double intersections. The grey boxes denote 
possible knottings and interlinkings with other portions of the loop (including possible additional intersections and kinks). In 
order to unclutter the figures, we only denoted partial connectivities, in all cases the loops have to be closed by connecting the 
open strands, with no particular restrictions, i.e. the resulting petal could have knottings and interlinkings with the rest of the 
loop. 
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FIG. 11. Schematic depiction of the first identity used to determine a and /3. 




/ < W{L^) > \ f < W{Lj) >\_ ^ -1 2 U < W{Lj) > \ 

[ < W{Lw^) > )- ^^^P^ 1^ < W{Lw) >)-'''^^['^k)[ a p)[< W{Lw) > J ' ^^^^^ 

We are interested only in computing < W{L^) >. We are not interested in the value of < W(L\y^) > since as we 
explained before it involves double lines, which would require a separate discussion involving extra skein relations (see 
subsection F). We do not know the values of a and /?. However, since we are computing the invariant as an expectation 
value of a Wilson loop, we expect it to satisfy the same relations that Wilson loops satisfy: the Mandelstam identities. 
We will see that imposing the Mandelstam identities is enough to determine the coefficients a and /3 uniquely and 
therefore to find the skein relation. 

We therefore need to discuss the content of the Mandelstam identities. As we mentioned before, these identities 
relate the behavior of reroutings of the loops at intersections and involve non-trivial information about the connectivity 
of the loop. We will restrict the explicit discussion to double intersections, but it is immediate to see that suitable 
generalizations can be found for more complex intersections. For double intersections, one can always consider a planar 
diagram, and two possible routings exist, as shown in figure |l^. Each routing leads to an identity. The identities, 

< Wivf'^) > + < VF(r;^24) > ^ < W{v§)W{vl^) > (37) 

< Wivl^'') > + < Wivl^^') > = < W{v]')W{v!^) >, (38) 



involve different loops, as shown in the figure ([1C|). Moreover, they are valid for a completely arbitrary loop that has a 
double intersection. In order to use these identities to allow us to fix the arbitrary parameters in the exponentiation, 
we will consider — for simplicity — its expression for a particular set of loops. We will then have to check that the 
resulting invariant is consistent with the identities for all possible loops. The loops we wish to consider are the ones 
obtained by reconnecting the loops of figure (|o|) with direct strands, i.e. adding no knottings or interlinkings and 
are shown in figures |ll| and The Mandelstam identities for these loops are, 

< VF(?7i) > + < W{m) >=< W{tjs)W{tj4) > (39) 
< VF(73) > + < T4^(74) >=< W{ji)W{-f2) > . (40) 
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FIG. 12. Schematic depiction of the second identity used to determine a and l3. 



We now make use of the skein relations for intersections with kinks (^,^). In equation (^9|) this allows to replace 
r]2 by the unknot and 773 and 774 by the two component unlinked link. Therefore < W{r]2) >= 1. 

In equation ( ^0| ) the use of the skein relations is shown in figure ^ to transform 73 and 74 to the unknot with a 
twist added. Using the skein relations ( |2^ , |2^ ) one gets, < 1^(75) >= and < 1^(76) >= 9'^. 

We now use (^6|) in combination with these expressions to determine the values of a and f3. It turns out we only 
need to use the expression of the exponential that appears in (js^) expanded up to second order in 1/fc. Using this 
expression to determine < W{r]i) >, < 14^(71)1^(72) > and < W{ri3)W{r]4) > one gets that a = 1 and /3 = 0. This 
completely characterizes the skein relation for the intersections (p6h. From there we conclude that, 



< W{Lw^) > 



g±4 (q=F4 — q^i) 
q^i 



< W{Li) > 

< W{Lw) > 



(41) 



From this expression, adding an subtracting with appropriate weights the expressions with the different signs, we 
can work out the usual skein relation for the Kauffman bracket without intersections. 



g3 < W{L+) > -q-- < WiL-) >= {q^ - q-^) < W{Lo) >, 
and also the definition of the expectation value for a planar straight through intersection, 

< W{Li) >= ^ , ^ _. [< WiL+) > + < W{L^) >]. 

(q4 + g 4) 



(42) 



(43) 



We therefore conclude that the expectation value of the Wilson loop in an SU (2) Chern-Simons theory is, up to a 
factor of 2 with the conventions of this paper, identical to the Kauffman bracket knot polynomial. 



C. Consistency for all loops and to all orders in 1/fc 

We have just shown that considering the exponential of the infinitesimal deformation to second order in 1/fc 
and requiring consistency with the Mandelstam identity for a particular set of loops uniquely fixes the values of 
the indeterminate coefficients a and f3 of the infinitesimal deformation. We need to check that the construction is 
consistent to all orders in 1/fc and for all possible sets of loops with planar intersections. We will now show that this 
is the case. In order to do this, the finite version of the form of the skein relation (Q), which again can be derived in 
the same form as the finite skein relation we derived above is useful, 

< W{L±) q^i < W{Lo) > -q"^i < W{Lq) > . (44) 

As we mentioned before, this relation depends on the connectivity of the loop, we are assuming it is such that Lj 
has one independent component, as shown in loop 77)^'^^ of figure (|To|). Notice that by subtracting the + and — sign 
versions of (Q) we obtain (|4^). The sum contains the information needed to derive the Mandelstam identity from 
the skein relations. Combining ( ^ ) with (^^ and (^,^) we get, 

< W{Li) > + < W{Lw) >=< W{Lw) > ■ (45) 
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/ w 

FIG. 13. The non-planar intersections we consider. The top row are three dimensional views and the bottom row are the 
associated crossings that arise in the skein relations. 



Notice that we are stretching the notation here, since this identity is not purely local, a given connectivity of the 
loop was assumed to derive it, as mentioned above. Therefore one has to include the connectivity we mentioned above 



in its interpretation. With that connectivity, equation (Hq) is exactly the Mandelstam identity 



the proof. One could have chosen a different connectivity when deriving 
Mandelstam identity (38). 



|) and then one would have arrived at the 



This completes 



D. Non-planar intersections 



Traditionally, in knot theory invariants have been formulated through planar projections of knots. When general- 
izations to intersecting knots are considered one may need to consider non-planar intersections as inequivalent. This 
has already been noticed for triple intersections [ p5[ . In this section we will discuss some skein relations satisfied 
by non-planar intersections. There are many types of non-planar intersections, we will restrict ourselves to some 
examples that have three of the four strands at the intersection in a single plane. We will show that there exists 
a one-parameter family of regularizations of the expectation value of the Wilson loop that is compatible with the 
Mandelstam identities for the case of non-planar double intersections. For one particular value of the parameter the 
intersections behave as planar ones. For other values of the parameter the intersections acquire common elements 
with under- and over-crossings and therefore imply the existence of distinct nontrivial generalizations of the usual 
invariants for the case on non- planar intersections. 



We consider a non-planar intersection as shown in figure 13, which we represent with two new types of crossings, 
labelled Lj and with obvious counterparts with a minus sign. 

In order to compute the skein relations we consider a deformation of a straight-through planar intersection as 
shown in figure |l^. The resulting integral is exactly the same as equation (|28|). The difference comes from the vector 
u. That vector vanishes in the z > Q part of the loop. This implies that when one wants to rederive equation 
( p9| ) one encounters ambiguities of the type S{y)Q{y) exactly as when we deformed an Lw intersection in a direction 
perpendicular to the plane. The way to handle this is again to introduce an indeterminate parameter a. The resulting 
expression then reads, intersections with 

ITT ^XIT 

dtSsit) < Tr(i7i2(7)i?34(7)) >= -a[-TT{H,2h)Hsi{l)) - — Tr(i7i2(7))Tr(i734(7))] (46) 







this can again be reinterpreted as, 

J dt6s{t) < W{Li) >= a[— < W(Lw) > < W(Li) >]. (47) 
In order to exponentiate, since the action mixes Lj and Lw we need to compute, 

( < W{Lf) > \ r ( < W{Li) >\_^^^ (^^\ (-a2a\( < W{L,) > 
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FIG. 14. The deformation of a straight through planar intersection into a crossing of the type 



and one can show that in order to have consistency with the Mandelstam identity, /3i = /?2 = a. Notice that a is 
undetermined by the Mandelstam identities and we therefore have a one-parameter family of definitions of the skein 
relation for the crossing. Exponentiating explicitly and identifying the variable k as before, we get, 

< W{Lf) > = g=Ff < W{Li) > - g=Ff ) < W{Lo) > (49) 

< WiL^) > = < w{Lo) > . (50) 

If a = we get that Lf = Lj and = Lw, so the non-planar and planar intersections are treated in the same 
way. 

It is also immediate to prove from the skein relations ( ^9|j4^ , ^ ) that, 

< W{Lf) > - < WiL+) > (51) 

< WiLf) > = < WiLi) > (52) 
< WiLf) > = < W{L-) > (53) 

so we see that for different values of the free parameter associated with the non-planar intersections we can have them 
play the same role as over-crossings, intersections and under-crossings. This highlights the relation between the value 
of the a parameter and the different types of regularizations it implies for the intersections. 

As in the non-intersecting case, all the regular invariant information of the < Vl^(7) > is concentrated in a multi- 
plicative "phase factor". One can divide by it and construct an ambient isotopic invariant, which for SU{2) is the 
Jones polynomial. In the non-intersecting case, the writhe ^(7) can be computed by evaluating the expectation value 
of a Wilson loop in a U{1) Chern-Simons theory, 

< W{^) > 1^(1) = g-(7), (54) 

with skein relations, 

w{L±) ^ ±1 (55) 

w{L±) = ±1 (56) 

w{Li) = wiLo) = (57) 

w{Lf) = ±a (58) 

If one now defines a polynomial J(7) through 

J(7) = ig-t-(7) < wi^) >, (59) 

the result is ambient isotopic invariant. 

The above definition of the Jones polynomial is also valid for multiloops, with a suitable generalization of the 
Abelian calculation of the writhe for multiloops. From the expression of the resulting Jones polynomial one can get 
an expression for the Gauss linking number of two loops with intersections. The resulting expression is consistent 
with lattice definitions of the linking number with intersections Pq,E7|. 
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E. Triple and higher intersections 



We will not discuss in detail the generalization to triple intersections in this paper, we sketch in this subsection how 
does one perform the generalization of the construction to that case. In the case of triple intersections, there are many 
types of possible independent vertices. In the case of planar intersections it is easy to see that one can deform them to 
double intersections, using the same techniques as for the Lw^s. For the case of non-planar triple intersections, there 
are 10 independent vertices. All of them can be related to double intersections through deformations similar to the 
one that connects Lj with L±. In order to exponentiate the infinitesimal deformations, one again has to exponentiate 
a matrix that connects the different intersection types. It will be a sparse 10 x 10 matrix. Again regularization issues 
will leave many coefficients undetermined and one would restrict them using the Mandelstam identities. It is not clear 
if the resulting polynomial will be completely determined by the Mandelstam identities or if new free parameters will 
appear. As in the double case, the Mandelstam identities are non-local and there are now three different possible 
connectivities of intersections that are needed to implement the identities. 



F. Loops with multiple lines 



Throughout this paper we have assumed we were dealing with loops that have each line traversed only one. If 
one would like the polynomial that is being derived to take values on the complete set of loops that is of interest in 
quantum gravity and gauge theories, one needs to consider the case of loops with multiply traversed section. This is 
also of importance if one is to view the integral of the exponential of the Chern-Simons form as a rigorous measure on 
the space of connections modulo gauge transformations. We present here a brief discussion of how could one consider 
double lines, but a complete description again requires further study. 

Let us consider a loop 7 and study the Wilson loop along 7^ =707. Consider a generic direction in space w 
such that the tangent vector to 7 is never parallel to w. Consider an continuous infinitesimal displacement of all the 
points of 7 along w such that an arbitrary point on 7 is kept fixed (we will call this point the origin of the loop). 
This produces a second copy of 7. If one deforms back and forth along w in this way one ends u p w ith two copies 



of 7 that are connected at the origin through an intersection of the type Lj, as shown in figure 15, Because this 



is a planar deformation, the value of the expectation value of the Wilson loop does not change. We have therefore 
reduced the problem of computing the expectation value for a loop traversed twice to the problem of computing 
the expectation value along a simply traversed loop with an intersection of one of the types studied before. It is 
worthwhile noticing that applying the Mandelstam identity (37) at the resulting intersection one gets the identity 
< W{j)W{j) >=< W{j o 7) > +2. 

If one wants to consider a more general situation, in which a loop could have portions that are traversed twice, 
one can extend the above result in a relatively straightforward manner. Even if there are planar intersections of the 
multiply traversed segment with another segment one can separate them without additional contributions. However, 
if there are triple non-planar intersections involving sections multiply traversed, the deformation will give non-trivial 
contributions. These contributions can be evaluated explicitly given a specific loop. 



V. CONCLUSIONS 



We have shown how to use variational techniques to obtain exact expressions for the knot invariant associated with 
the expectation value of the Wilson loop in Chern-Simons theory. The method is completely general in the sense that 
it can be used in Chern-Simons theory with any semisimple gauge group with small generalizations. We have worked 
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out explicitly the value of the invariant not only for smooth loops but also for loops with double planar and non-planar 
intersections. The resulting invariant is compatible with the Mandelstam identities of the gauge group and therefore is 
suitable for providing invariants of interest as quantum states of topological field theories and quantum gravity. Having 
a well defined linear function on the space of loops compatible with the Mandelstam identities may allow also, using 
the techniques of Ashtekar and collaborators to define in a rigorous way a measure in the space of connections 
modulo gauge transformations dii(A). Such a measure would allow to give rigorous meaning in a mathematical sense 
to expression of the form / dAe''^^''^ f{A) for any gauge invariant function f{A). The generalization of the work 
described in this paper to triple intersections is straightforward and relevant for quantum gravity applications. The 
possibility of computing explicitly the knot polynomials associated with Chern-Simons theory for any group is clearly 
of relevance in other physics applications, like the recent discovery of Chern-Simons states in supergravity has shown 
p8[ . It is expected that these results will be extendible to A'' > 1 supergravity where this method will provide with a 
simple method of characterizing the potentially new invariants that may arise. 
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APPENDIX A: COMMUTATIVITY OF THE 5 OPERATORS 



We wish to evaluate the successive action of two 5 operators on < W{'j) >. We only discuss the case of a regular 
(non-intersecting) point of the loop. The first 6 acts as indicated in (19), and using the Fierz identity we get, 

Suit) < W{j) >=-^ f dsr,-{s, tW\s, t)eabc / dy-6^{f^{s, t) ^ y) < W{^) > . (Al) 
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The second S has two contributions, stemming from the action of the loop derivative on the loop dependence of 
< Wij) > and respectively. The term that contributes to the commutator is the latter, since the contribution on 
the Wilson loop is the same in both orders. Let us therefore evaluate explicitly the non-trivial contribution, 

[5{ti),S{t2)] ^ — — / dsi / ds2f]''isi,ti)rf isi,ti)ri''{s2,t2)r]' is2,t2)eab[cdd]S^ix - rj{si,ti))\^^r,{s2A2) 

-{h^t2) (A2) 

= r- / dsi / ds2f]°-{si,tx)r]' {si,ti)r]''{s2,t2)eabcds2S^{il{s2,t2) - visi^h)) 

^ Jo Jo 



+ — J dsi J ds2'n''{si,ti)r/ {si,ti)r]' {s2,t2)eabddt2S {v{s2,t2) - v{si,ti)) 

-{h^t2). (A3) 

One can now pick a coordinate chart in which S^{'r]{s2,t2) — rj(si,ti)) = Z5{s2 — si)S{t2 — ii), where Z is a 
regularization dependent factor that does not depend on S2 or ^2- Using the distributional identity, 

f{y)dy5{y - x) = f{x)dy6{y - x) - 5{y - x)d,f{x) (A4) 

it is immediate to check that the contributions vanish. The result is obviously dependent on a regularization choice. 
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